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Schemes for topological quantum computation are usually based on the assumption that the 
system is initially prepared in a specific state. In practice, this state preparation is expected to 
be challenging as it involves non-topological operations which heavily depend on the experimental 
realization and are not naturally robust against noise. Here we show that this assumption can 
be relaxed by using composite anyons: starting from an unknown state with reasonable physical 
properties, it is possible to efficiently distill suitable initial states for computation purely by braiding, 
i.e., reversible gates. This is achieved by encoding logical information in a subsystem code with gauge 
system corresponding to the internal degrees of freedom of composite objects. 



I. INTRODUCTION 

The effects of imperfections and decoherence are ar- 
guably the most serious obstacle faced when trying to 
design a quantum computer. Their study has led to the 
development of elaborate schemes for fault-tolerant quan- 
tum computation. These efforts have culminated in the 
threshold theorem [l|, [24| |32| , which establishes that re- 
liable computation is in principle possible given a suffi- 
ciently low error probability of each individual gate. A 
central idea underlying these results is the concatena- 
tion of codes, first formalized by Forney [l5[ and later 
generalized to quantum error-correcting codes [ltH fl7j . 
Here, logical information is encoded in a recursive fash- 
ion, where the qubits constituting a codeword are them- 
selves represented by an error-correcting code. This leads 
to an exponential reduction of the effective error rate in 
the number of concatenation levels. 

Topological quantum computation (TQC) is an alter- 
native (and in some sense complementary) approach to 
fault-tolerance [H H HH, H • The basic steps for TQC 
can roughly be summarized as follows (but see Q for an 
alternative complete set of basic operations) : (i) create a 
certain initial state of non-abelian anyons (e.g., the state 
of many particle-antiparticle pairs, as explained below), 
(ii) braid the anyons around each other to perform gates, 
and (iii) measure the resulting state (that is, the total 
topological charge of a region). The non-local nature of 
the encoding used here is expected to lead to significantly 
lower error rates compared to more standard implemen- 
tations. Moreover, gates are intrinsically robust due to 
the fact that (ii) only involves topological operations. 

Unfortunately, steps (i) and (iii) of TQC are gen- 
erally system-dependent, and require performing non- 
topological operations which are not naturally protected 
against noise. For example, measurement of topological 
charge involves making interferences from the observa- 
tion of a non-topological quantity (e.g., current) through 
approximate relations valid only in the given experimen- 
tal setup. Here we focus on the state preparation prob- 
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lem (i). We give a scheme for creating a suitable ini- 
tial state purely by braiding, starting from an unknown 
state supported on a physically relevant subspace. This 
reduces TQC to a sequence of braids followed by mea- 
surement. It eliminates all unprotected operations apart 
from those at the end of the computation. 

The basis for our scheme is the fact that a region of 
a topologically ordered 2D-quantum medium can be re- 
garded as a dot where a composite anyon resides. This 
composite object carries the total (topological) charge 
of all the anyons contained in that region. Anyon sys- 
tems therefore provide a natural mechanism for concate- 
nated coding: fusion basis states of composite objects 
can be used to encode information. At each level of 
concatenation, a number of subregions are joined into 
a region forming a new composite object. Ignoring the 
degrees of freedom within subregions, the resul ting en- 
coding has the structure of a subsystem code @, I22l l23j. 
where the Hilbert space decomposes as % = (C <g> V) © 
and the logical information is encoded in C. Here T> de- 
scribes the internal degrees of freedom of the composite 
anyons. Importantly, computations in the code space C 
can be performed in the same manner as for the bare 
anyons. In particular, the complexity of braiding two 
regions around each other only scales polynomially with 
the size of their boundary. (This can be seen explic- 
itly for schemes based on local code deformations.) A 
similar statement can be made for the measurement of 
composite topological charge, which can be achieved by 
anyon interferometry [6|, for example. To compute, it is 
therefore sufficient to be able to prepare any state sup- 
ported on the subspace \^i n it) <8> £>, where \^mit) is a 
state of a fixed number of anyons representing the logical 
state |0)® n of n computational qubits. 

We consider (mixed) states supported on a physically 
relevant subspace. We show that any such state p can 
be turned into a state supported on \^i n u) <8> T> using 
reversible gates. Moreover, we construct braid sequences 
implementing this efficiently. In other words, we con- 
struct a distillation procedure which concentrates the en- 
tropy of p into T>. Using this distillation procedure cir- 
cumvents the need for measurements when preparing the 
initial state for computation. 

The remainder of this paper is structured as follows: 
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In Section UH we give some background on anyons and 
explain the action of the braid group. In Section |III[ we 
present our scheme for distillation of composite anyons. 
We conclude in Section HVl 



II. ANYONS AND THE BRAID GROUP 
ACTION 

A. FQHE and qudit lattice Hamiltonians 

The fractional quantum Hall (FQH) effect is the most 
well-known example of a topologically ordered state of 
matter. At filling fraction v = 5/2, this system is be- 
lieved to be in the universality class of the Moore-Read 
Pfaffian state [30| . which is the k = 2 state in the 
Read-Rezayi-sequence of states [13, Hl|. Similarly, for 
v = 12/5, the system is believed to be described by 
the k = 3 Read-Rezayi-state. Quasi-particle excitations 
above these states are superconducting vortices carry- 
ing fractional (electric) charge e/4. They are predicted 
theoretically to correspond to the anyons described by 
su(2)fc Chern-Simons-Witten theory. In principle, these 
realize two of the most basic non-abelian phases: the 
Ising TQFT (for k = 2) and the Fibonacci TQFT (for 
k = 3). At present, though, there is only some experi- 
mental evidence for the existence of these quasiparticles, 
and it is restricted to the former case [lCj, [3J, |42j . 

FQH systems are not the only candidate systems for 
building a topological quantum computer. Of particular 
interest are qudit lattice Hamiltonians with topological 
order. Examples are Kitaev's toric code Kitaev's 
honeycomb model [2(| or Levin- and Wen's Hamilto- 
nian [27[ . While these Hamiltonians are unlikely to occur 
in nature, it is in principle possible to engineer systems 
with corresponding interactions, and some experimental 
proposals for doing so exist (see e.g., [HI, [29j]). With an 
appropriate choice of Hamiltonian, anyons corresponding 
to any "doubled" (non-chiral) TQFT can be obtained 
from such a construction. This approach to achieving 
topological order also has an additional advantage: the 
local Hamiltonian terms can be interpreted as local stabi- 
lizers of a quantum error-correcting code of qudits. Using 
this code for storage and processing of quantum compu- 
tation merely requires the ability to measure syndromes 
locally, instead of requiring a system described by the 
Hamiltonian. In this setting, the problem of decoher- 
ence at non-zero temperature (see e.g., 0) can be over- 
come by performing active error-correction Q. Indeed, 
this gives rise to excellent threshold results for quantum 
error-correction [35l l36j . 

To use a topologically ordered medium for computa- 
tion, one requires the experimental ability to manipu- 
late, i.e., create, move (braid) and measure anyons. For 
FQH states, proposals for implementing anyonic charge 
measurements are most developed. These employ beams 
of probe quasiparticles (e.g., constituting a current along 
an edge) , which are sent through an interferometric setup 



(e.g., by deforming the edge around a region to be mea- 
sured). Because of the non-abelian statistics, topological 
charge information can be deduced from the interference 
pattern observed by detecting whether a quasiparticle 
is present at the outputs of the interferometer. Since 
the latter detection requires coupling to e.g., the electric 
charge of the probe anyons, this is a rather indirect pro- 
cedure which is susceptible to noise. Nevertheless, this 
is at present the most promising experimental approach 
to observing anyonic statistics. In contrast to measure- 
ments, current theoretical proposals in the FQH setting 
for other operations, such as controlled movement of in- 
dividual particles, perhaps have little chance of actually 
being realized. 

The manipulations required to perform a topological 
quantum computation in a qudit lattice system consist of 
local unitaries and measurements acting on a small num- 
ber of qudits at a time. This has been studied in detail for 
the toric code [l{| and the Levin- Wen- Hamiltonians [27| ■ 
Anyons are created and moved by applying certain "rib- 
bon operators" (which can be decomposed into locally 
acting operators). Some proposed schemes (l3l [25| di- 
rectly use a topologically degenerate ground space and 
proceed by adiabatically deforming the Hamiltonian in 
a local fashion. Measurements of topological charge re- 
quire decoding the encoded information, and thus are 
automatically subject to noise. Corresponding circuits 
are discussed e.g., in [9l [l2l |2(|. 

B. Fibonacci anyons 

We formulate our scheme for the Fibonacci theory, the 
simplest computationally universal anyon model. We 
proceed with a short description of the corresponding 
braid group representation. For more details, we refer to 
the literature. Excellent introductions can be found e.g., 
in 

The Fibonacci model has only one non-trivial parti- 
cle r with fusion rule r x r = 1 + t. For simplicity, we 
consider a (infinite) plane with N localized dots, each of 
which can support a r particle (but may also have trivial 
topological charge). We assume that the total topologi- 
cal charge is 1. The Hilbert space then decomposes into 

N ( N ) 

a direct sum H = ©„ =0 ©i=i Ki> where V£ is the space 
of n r-anyons localized at dots specified by the subset i. 
Universal quantum computation can be performed in a 
subspace C N ^ 4 C =: Vn, but here we take a more 
general viewpoint. 

The space of N r-anyons is characterized by the 
action of the braid group. In the FQH setting, it is the 
degenerate eigenspace of states with localized excitations 
at specified locations (i.e., the dots). In topological com- 
putation schemes based on the degenerate ground space 
of a qudit lattice model, the space is a certain sub- 
space of the ground space. Here dots correspond to de- 
fects, i.e., holes in the latttice structure or (possibly ex- 
tended) areas where the local Hamiltonian terms are dif- 



ferent from the bulk. The mapping between states of 
the physical system and the abstract anyonic space Vk 
is discussed extensively in the literature (see e.g., [33l.l37| 
and [l9|,[2(J[2l[ for lattice systems). Similarly, braids cor- 
respond to certain (physical) operations on the system, 
as discussed above. Here we will restrict our discussion to 
the abstract anyonic formalism which provides a system- 
independent description of the action of the braid group. 

We use the standard diagrammatic representation of 
states and operators on H (see e.g., [|| for more details): 
To every tree with trivalent vertices and N leaves, we 
can associate an orthonormal basis of % with elements 
given by different labelings of the edges by {l,r} consis- 
tent with the fusion rules (i.e., with the number of edges 
labeled t incident to any vertex not equal to 1). We fix 
the label at the root of the tree to 1 corresponding to 
the trivial topological charge at the co-boundary. The 
labels attached to the edges connected to the leaves de- 
termine which of the spaces a labeling corresponds to. 
The tree specifies the order of fusing, and basis elements 
corresponding to different trees are related by so-called 
F-moves. In the following diagrams, we represent edges 
with the label r by solid lines, whereas edges with the 
trivial label 1 will be dotted or omitted altogether. An 
F-move then is a local substitution of subgraph by a lin- 
ear combination as specified by the identities 
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where <p — v ^ +1 is the golden ratio. 

The braid group Bn on N strands is generated by 
the elementary braids <t\,... , crjv-i satisfying the Artin 
relations (Ji(Jj = <Jj(Ji for \i — j\ > 2 and <Ji(Ti + i<Ji — 
ai+i(Ti<Ji + i for I < i < n — 2. To describe the action 
of Bn on %, one represents the generators as 

) 



ii + 1 N 



ii + 1 N 

Group multiplication then corresponds to stacking such 
diagrams. To apply a braid b £ Bn to a basis ele- 
ment |^) G TL, one attaches the graphical representation 
of b to the leaves of the tree specifying |^}, and then 
reexpresses the result in the original basis using a set of 
local rules. These allow to arbitrarily add and remove 
edges with the trivial label 1 in addition to isotopic de- 
formations that keep the leaves fixed, F- moves as in (pj, 
substitutions of bubbles by scalars as specified by 
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FIG. 1: Left: A possible arrangement of the dots (small cir- 
cles) in the plane; we write basis states arranged from left 
to right, but the actual dots need not be arranged on a 
line. A typical state in H p hys has some TT-pairs (depicted by 
filled circles), and trivial composite topological charge for any 
larger region containing a pair of dots (dotted lines). Right: 
The result of applying the described braiding transformation 
(wo&)® 6 to each group of four is a state with composite anyon 
pairs in the indicated locations. 



and resolution of crossings using the i?-matrix. The lat- 



ter take the form 
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The space V£ has dimension dimV^* = 
\4' n ^ 1 — (— ■ It is invariant under the 

action the purebraid group PBn on N strands. This is 
the subgroup of the braid group Bjy where each strand 
begins and ends in the same position. The restriction of 
this group action to V£ has a dense image in the unitary 
group PU (dim V£) modulo global phase factors [14| 
(i.e., is computationally universal). More precisely, 
PU(M) = SU(M)/Z M , where Z M is the subgroup of 
SU(M) generated by e 27 "/ M • id. 



III. DISTILLING COMPOSITE ANYONS 
A. The setting 

Let us identify physically relevant states. We envision 
the dots to be arranged in spatially separated pairs, see 
Fig. [1] Let \vac) € Vq be the (up to a phase) unique 
state of no particles. We could imagine the system to be 
initially prepared in this vacuum state, but subsequently 
undergo processes which create particle-antiparticle pairs 
locally at pairs of dots. We expect tunneling of (topolog- 
ical) charges between different pairs to be exponentially 
suppressed in the inter-pair distance (cf. [j|). Due to the 
geometric arrangement of the dots, the resulting state 
will therefore essentially be a mixture of superpositions 
with trivial charge for each pair of dots, but having tt- 
pairs within some of them. 

It is natural to ask for the distribution of rr-pairs in 
such a state. Here we will make the minimal assumption 
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that within region Q containing mu pairs of dots, there is 
a non-zero number of rr-pairs except with a probability 
exponentially small in mn- More precisely, let Pq de- 
note the projection onto the subspace spanned by states 
where all pairs in ft have trivial total topological charge 
and at least one rr-pair is present. We assume that the 
states are supported on the image of Pq up to a fraction 
exponentially small in ran. States of this form are the 
starting point for our braiding procedure. 

As a prototypical example of a state satisfying our as- 
sumptions, suppose p = $p N ^ 2 {\vac)(vac\) is the result 
of a process $ v which creates particle pairs with proba- 
bility p locally. Then the number of pairs in a region Q 
is binomially distributed with parameters (mn,p), and 
tr(P n p) = l-(l-p) m ". 

We will assume that after the initial pair creation pro- 
cesses, the dots are moved apart from each other so 
that further particle tunneling is suppressed. We will 
show how states satisfying our assumptions can be effi- 
ciently mapped to a subspace \^i n it) €5 D (up to a neg- 
ligible error) by braiding, where \^i n it) £ Vik is the 
state of k composite rr-pairs each with trivial topological 
charge. The latter can be identified with the state [0)® fe / 2 
of k/2 logical qubits in a standard encoding. 

We conclude this section by commenting on the phys- 
ical systems where our scheme may be applicable. Be- 
cause of the restricted repertoire of currently realizable 
operations in the FQHE setting, the proposed procedure 
may be more useful in the lattice Hamiltonian setting (if 
such a system can be realized). Furthermore, in quantum 
Hall states, anyons typically carry additional quantum 
numbers such as electric charge. This should simplify 
their detection by local measurements; it also leads to 
an additional energetic cost for having anyon pairs. In 
contrast, anyons encoded in the ground space of a lattice 
system (as in [IHHBj]) have degenerate energies even for 
different topological charges. Also, these anyons are not 
naturally localized, but correspond to topologically non- 
trivial loops around holes. To measure these topological 
charges requires non-local measurements which are sus- 
ceptible to noise. Thus such lattice systems are natural 
candidates for the application of the proposed scheme. 



B. Distillation of a single pair 

For simplicity, we first consider the problem of "dis- 
tilling" a single composite rr-pair, i.e., k = 1, and gen- 
eralize to k > 1 later. We show the following: Con- 
sider a region SI = £l A U Ci B partitioned into disjoint 
sets £l A ,Q B each containing toq/2 = ttiqa — itiqb pairs 
of dots. Then there is an efficiently computable se- 
quence of poly(mn, log 1/e) elementary braiding opera- 
tions within il with the following effect: the image \ty') 
of any state \&) G PqH. can be approximated in trace 
norm to distance e by a superposition of states with a 



composite rr-pair between H, A and f2 B . Pictorially, 




where the notation P» e indicates that the result is ap- 
proximate up to an error s in trace norm. This pic- 
ture is to be understood in a basis where fl A and Sl B 
are the leaves of two subtrees joined together. In other 
words, the resulting state is supported on the subspace 
l^init) ® 2? c C ® V, where the code space C = V% = C is 
the space of 2 composite r-anyons defined by the subre- 
gions f2 A ,0 B , and the space V corresponds to the inter- 
nal degrees of freedom in the subregions (represented by 
boxes). Elementary braiding operations are defined as 
clockwise or counterclockwise exchanges of regions con- 
tained within f2, and thus have a polynomial complexity 
in the length of the boundary of fl. 

To establish our main claim, consider the special case 
where vtlq = 2, i.e., fl A and Q B each contain a pair of 
dots. Here the subspace Pq H is spanned by states which 
locally have the form of the three states 

V y,y y,y y. 

Our goal is to map such states to states of the form 

1*1 (4) 

by an appropriate braid of the dots. By the density result 
mentioned above (applied to the space V4), there is an 
element b S PB4 which approximately maps 

y y. • ^ (5) 

up to an irrelevant phase. Using the fact that 6 is a pure- 
braid, we conclude that the other two states are mapped 
to themselves up to phase factors, i.e., 

V v >'•■' y y, y • < ' \j y 

(6) 

Note that, by the Kitaev-Solovay theorem 2lj], a se- 
quence of C(log l/e) a (where C is a constant and a « 4) 
generators implementing a braid b with the desired prop- 
erty can be found efficiently. Executing this sequence is 
the first step in our procedure. It remains to map the 
two states in © to states of the form Q, while leaving 
the rhs. of §5§ invariant. For this purpose, we use the 
injection weave construction proposed in Q. This is an 
element w G -B3 with the following properties: (i) up to 
an error e, it acts as the identity on the space of three 
r-anyons, (ii) it is a weave, i.e., only a single strand, 
the warp-strand is moved around 2 stationary strands 
and (iii) the warp strand starts in the right-most posi- 
tion, and is transferred to the leftmost position. Fig. [5] 
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FIG. 2: The form of an injection weave. The warp strand is 
shown in bold and the box stands for a sequence of braids 
whose details are irrelevant here (explicit constructions of 
such weaves are given in ||, HU). The injection weave acts 
as the identity on the space of three r-anyons up to a small 
error. 



shows the form, i.e., conditions (ii) and (hi) defining 
such a braid. As for the braid 6, an injection weave w 
of length C(logl/£) Q can be found efficiently using the 
Kitaev-Solovay algorithm (see [39| for a general discus- 
sion of the existence of injection weaves and its appli- 
cation to the construction of "pureweaves" for universal 
computation). It is clear that applying such an element w 
to the three rightmost strands maps 
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for some phases 4>\,4>2, while approximately preserving 
the state on the rhs. of ([3]) up to an error s by prop- 
erty (i). Importantly, the precision e of b and w can be 
improved exponentially by a polynomial increase in the 
length of the braids. We conclude that the composition 
wob achieves the desired result in the case where mo = 2. 



C. Distillation of multiple pairs 

To prove our claim in the general case, consider a re- 
gion O containing = 2 pairs of dots. We divide 
these into groups of 4 dots (two pairs), and apply wob to 
each of these groups. Picking an area around each pair, 
we can think of the result in terms of "superdots" (each 
containing a pair) where composite r-anyons may reside. 
Inspecting the previous analysis, we conclude that by go- 
ing from dots to composite objects, we half the number 
of objects, while effectively increasing the density of rr- 

be defined in the same 



pairs. More precisely, let P n 
way as Pq , but for the 2 
map (w o 6)® 2 takes 



1 pairs of superdots. Then the 

|*') e P^n , (7) 



I*) 6 PnU ^2e-2'- 

and is composed of 2C(log 1/e)" ■ 2 l ~ x elementary braid- 
ing operations. Here we used the triangle inequality to 
give an upper bound on the error. 

Finally, we use this procedure recursively: at the r-th 
level (with r decreasing from £ to 1), we divide the given 

2 r+1 regions into groups of four, apply (w o b)® 2T 1 and 
obtain 2 r superregions containing pairs of regions. This 
is shown schematically in Fig. [3] (cf. Fig. [l| . We are inter- 
ested in composite rr-pairs among pairs of these super- 
regions (at most 2 r_1 may be present), and define -Pq -1 



w 



FIG. 3: The first two levels of the composite braid used for 
the preparation of the initial state. Thick strands correspond 
to composite anyons. 



accordingly. By induction using ([7]), we conclude that 
after I — 1 steps, we obtain a map 

|*) EPnn >=h5 |*") e PqH with precision 

i 

5={2e)Y,y 1 <2^(2e) . 

r=l 

Furthermore, this requires 2C(logl/e) Q ■ Ylt=i 2 r_1 < 
2 •2C(log l/e) a elementary braiding operations. We con- 
clude that any desired precision s can be obtained with 
run- poly (log ma, log 1/e) elementary braiding operations. 
Furthermore, by definition of P^, the resulting state has 
a (composite) rr-anyon between two regions of Q hav- 
ing 2 l points each. We have thus established our main 
claim for k — 1 encoded rr-pair. 

Let us now generalize this result to the generation of 
k > I composite rr-pairs. Since the transformation Q 
does not affect the region outside ft, we can simply re- 
peat it k times. That is, we fix a partition Uq=i ^« of 
the system into sufficiently large regions £l a , each con- 
taining m pairs of dots (m is a power of 2). As argued 

above, the space H p h y s '■= (IIa=i -^fO^i which is the 
space of states having at least one rr-pair in each region, 
essentially supports all states of interest. Applying the 
transformation ([3]) for each a maps any such state into 
a state of k encoded rr-pairs with error at most ke, i.e., 
(for k = 3) 

The code space C = Vik is the space of 2k composite r- 
anyons defined by subregions Clf, Slf , . . . , f2^, 0^, where 
we have chosen a bipartition f2 Q = U for each 
a = 1, . . . , k with rriQA — m^B = m/2. The space T> rep- 
resents the internal degrees of freedom in the subregions. 
Note that the resulting state can be used to represent the 
state |0}® fc / 2 of k/2 logical qubits in a "computational" 
subspace C fe / 2 C V2k using a standard encoding (28j . 

In summary, our braiding sequence is efficiently com- 
putable and creates k composite rr-pairs from any 
state supported on Hphys — (Ila=i )?^- Any de- 
sired precision e is obtained using a polynomial num- 
ber poly(m, fc, log 1/e) of elementary braiding operations. 
The latter are defined as clockwise or counterclockwise 
exchanges of regions contained within the sets fl a , and 
thus have a polynomial complexity in the length of the 
boundary of fl a . 
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IV. CONCLUSIONS 

Wc have demonstrated that using composite anyons in 
conjunction with specific braiding sequences can elim- 
inate the need for preparing a particular initial state 
for topological quantum computation. Our scheme is 
highly parallelizable, and bears a remarkable similarity 
with entanglement renormalization methods (4lj . sug- 
gesting that these concepts might extend to the set- 
ting of anyons. Our work is merely a proof of principle 
and is restricted to Fibonacci anyons, though our argu- 
ments can be adapted to give procedures for the anyons 
of sw(2)3 Chern-Simons-Witten theory along the lines 
of [18|, as well as to e.g., the doubled Fibonacci the- 
ory. Whether or not similar schemes can be constructed 
for more general anyon models is an interesting open 
problem which goes beyond commonly considered uni- 
versality questions. Computational universality allows to 
implement NOT-type gates using pure-braids, and also 
gives rise to injection weaves (§1. l39|. but this is not nec- 



essarily sufficient to provide a distillation scheme along 
the lines explained here. Another challenging problem 
is the examination of the robustness of the scheme with 
respect to different initial states having composite topo- 
logical charge. In practice, it may be necessary to use a 
combination of the techniques presented here and more 
traditional methods relying on measurements and error- 
correction. 
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